Lateral organization in mixed lipid bilayers supported on a geometrically patterned 

substrate 
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The organization of lipids in biological membranes is essential for cellular functions such as signal 
transduction and membrane trafficking. A major challenge is how to control lateral lipid composition 
in supported membranes which are crucial for the design of biosensors and investigation of cellular 
processes. Here, we undertake the first theoretical study of lateral organization of mixed lipids in 
bilayers induced by a geometrically patterned substrate, and examine the physical mechanism of 
patterned substrate-induced structural formation in the supported lipid bilayers. A rich variety of 
composition segregations of lipids are regulated, and the results can account well for most recent 
experimental works [Yoon, T. Y. et al. Nat. Mater. 5, 281(2006) and Parthasarathy, R. et al. 
Langmuir 22, 5095(2006)]. The present study provides a comprehensive understanding of mechan- 
ically controlling the spatial organization of membrane components by unifying these experimental 
evidences. 
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Introduction 

Biomembranes are self-assembled bilayers of lipid 
molecules. Recently, supported membrane [l|, Q has at- 
tracted extensive interest, not only due to its importance 
in studying the properties and functions of biological 
membranes, but also due to its potential applications in 
design and fabrication of biological devices. [1, |3j H [1, 0| 
Supported bilayers can maintain the structural and dy- 
namic properties of free bilayer such as the lateral fluidity, 
share many similarities with the natural membranes, and 
therefore are widely used as cell-surface models that con- 
nect biological and artificial materials. To realize the 
bio- functional supported membranes or design novelly bi- 
ological devices, we usually need to manipulate the sup- 
ported membranes in nanometer or micrometer scale, [g 
For example, one can apply a lateral electric fleld to re- 
arrange the distribution of lipids in a supported mem- 
brane for some related cellular processes. [8] Additionally, 
one can also separate a supported membrane into well- 
defined membrane domains by introducing some barriers 
into it, and micropattern it by the help of approaches 
of photolithographic patterning, microprinting, or mi- 
crofluidic flow patterning. 0, Uo] Supported membrane 
micropatterning has attracted tremendous attention in 
recent years because of its potential applications in the 
investigation of immunobiology, drug discovery, design of 
biosensors, etc. [3] 

Most recently, Lee and co-authors [Til [T^. [Tst presented 
an experimental approach to manipulate the lipid segre- 
gation in a mixed supported membrane by use of the ge- 
ometrical property of substrate. They studied the lateral 
organization of lipids in a mixed membrane supported 
on a substrate with a groove, [III and found that, due 
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to different effective molecular shapes of mixed lipids, 
the two kinds of lipids spontaneously segregate, where 
the lipids with big spontaneous curvature are preferen- 
tially localized at the grooves. In a most recent exper- 
imental work, Yoon et al.[l^ studied the effect of to- 
pographical substrate on the lipid raft formation in the 
supported membrane. In their system, the topography of 
the substrate is either nanocorrugated or nanosmooth, 
and the membrane supported on such a substrate will 
be curved at certain predetermined positions to follow 
the substrate structure. Because of large bending rigid- 
ity of sphingolipid/cholesterol-rich Iq domains, there ex- 
ists a large free energy when they stay in the curved 
regions. Thus, the macroscopic rafts which are coars- 
ened by the nanorafts, can emerge in the nanosmooth 
regions of substrate, whereas there are only nanorafts 
in the nanocorrugated regions. On the other hand, 
Parthasarathy et al.[T5| proposed a curvature- mediated 
modulation of phase-separated structures in membranes, 
and they found that, beyond a critical curvature of mem- 
brane geometry, the cholesterol-rich lipid-ordered lo do- 
mains stayed in the small-curvature regions, while the 
cholesterol-poor lipid-disordered Id domains were prefer- 
entially localized in the large-curvature regions due to 
their different rigidities. Besides the lipid domains in 
supported membrane systems, it was also found that 
membrane curvature can induce phase separation in the 
mixed giant unilamellar vesicles (GUVs) due to the dif- 
ferent rigidities of different domains. fldL [TgI. Iv^ 

Further research will be needed to uncover and re- 
veal all the possible lipid structures and physical mech- 
anisms behind experimental evidence from lateral orga- 
nization in supported membranes. As is well-known, it 
is still difficult to systematically probe and visualize the 
laterally heterogeneous structures of lipids from tens to 
hundreds of nanometers in size with current experimen- 
tal technologies, [l^ There is therefore an urgent need 
to gain greater theoretical insight into the physical pic- 
ture of how membrane organization is governed. Here, 
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we systematically investigate the lateral organization in 
a mixed bilayer consisting of symmetric lamella-forming 
lipids and asymmetric cylinder-forming lipids supported 
on a geometrically patterned substrate, and explore phys- 
ical mechanism behind structural formation of lipids by 
using self-consistent field theory (SCFT). We expect the 
present study to offer a universal and thorough strategy 
towards efficiently controlling the lateral phase separa- 
tion of mixed lipids in bilayers by varying the substrate 
roughness, in addition to yielding insights into several 
experimental findings. 



Model and Methods 

We consider a system composed of rii lipid species A 
and 712 lipid species B dispersed in iis hydrophilic ho- 
mopolymer solvents, [l9| which is supported on a geo- 
metrically patterned substrate(see Fig. 1). The head- 
group volumes of lipid species A and B are Vhi and 
Vh2(vhi < Vh2), respectively, and their tails have the 
same length consisting of Nt segments of segment volume 
Po^""^-!!^ The solvent chain consists of Ng segments with 
segment volume p^^ . In the present problem, we choose 
Nt = Ns = N without loss of generality. The conforma- 
tion of the ath chain of lipid A, lipid B, and homopolymer 
solvent can be characterized by continuous space curves 
ra,i(s), ra_2(s), and rQ_s(s), respectively, and all of the 
chains are completely flexible. [2l| We assume that the 
system is translational invariant in y-direction, then 
the surface of patterned substrate is characterized by a 
periodical function S{x), which is only dependent on x. 
In our calculation, we choose S(x) as follows, 



^(^^ |i/I(l + sin(^)), if n 



if n < ^ n + i. 



1. 

(1) 

Here, Hq is the depth of the grooves, Lq is the period of 
the substrate, and n{= 0, 1, 2, • • • ) is the number of the 
periods in the patterned substrate. 

To characterize the system, we define the concentration 
operators of headgroups $?ii(r) and tails $ti(r) of i-type 
lipids(i = 1,2 denotes lipid species A and B, respectively) 
and solvents $s(r) as follows. 



$t,(r) = — V / dsS{r - ra.i(s)) 

N "° 

$s(r) = — V/ dsS{r - rc,,sis)) 



(2) 
(3) 

(4) 



where pht{r) = ^ J^l^i ^i^ - i"a.i(l)) is the dimension- 
less number density of headgroups of i-type lipids, and 
7i = VhiPo/N. 
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FIG. 1: A schematic of lateral organization in two-component 
lipid bilayer supported on a geometrically patterned sub- 
strate. The lipid species A and B have different effective 
shapes (cylindrical- and cone-like) due to different head-tail 
symmetries. The substrate is periodically patterned, and each 
period is composed of a flat step and an arc-like groove with a 
certain curvature. The thickness of the lipid bilayer on the flat 
steps is Lm, the depth of the grooves is Hq, and the substrate 
period is Lq. 



We assume that the system satisfies the incompressible 
condition $,,i(r) + $ti(r) + *h2(r) + ^t2(r) + *s(r) = 
<i>o(r), and $0(1") is given by, 22] 
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< z < S{x), 
S{x) < S{x) + t, 
z ^ S{x) + e, 

(5) 

where e is a sufficiently small quantity. 

The SCFT method is very powerful to inhomogeneous 
macromolecular systems (see, for example, reviews psl 
[24 . 25]) and has been applied to investigate the lipid bi- 
layer or lipid bulk phase behaviors in recent years. [2^ [13, 
m, m, jsO] Leermakers and co-authors used a lattice ver- 
sion of SCFT to study the lipid-bilayer system, [H, [13] 
where the conformation of lipid chains is described by 
random-walk and RIS (Rotational Isomeric State) statis- 
tics. However, in most of their works, they assumed that 
the lipid bilayer was laterally homogeneous for compu- 
tational convenience. On the contrary, Schick and co- 
authors adopted the off-lattice version of SCFT to study 
the lipid svstems.[28l. [29l. [30| where the amphiphilic na- 
ture and chain stretching of lipids are still captured. 
Compared to the lattice version of SCFT used by Leer- 
makers and co-workers, the off-lattice SCFT can account 
reasonably well for lateral inhomogeneity of the lipid bi- 
layers, although less molecular details are considered. 
Here, we emphasis on the description of laterally hetero- 
geneous structures and phase separation of lipids due to 
the interplay among the conformational entropy of lipids, 
effective molecular shapes of lipids, and substrate geome- 
try. Therefore, the molecular details of lipids are not im- 
portant, and it will be more reasonable to use the coarse- 
grained model in the framework of off-lattice SCFT [3^ . 

In the framework of self-consistent field theory (SCFT), 
we first transform a hamiltonian of the lipid system into a 
coarse-grained field theory description by employing the 
simple Gaussian model for lipid chains and treating their 
interactions through mean-field Flory-Huggins form. The 
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partition function of the system is given by 
1 



Z = 



ni!7i2!ns! 



X 5($o(r) - *hi(r) - $ti(r) - ^^aW - 4t2(r) 
- $,(r)) X exp{-po / dr[x{^hiir) + $h2(r) 



(6) 



+ $,(r))($ti(r) + $t2(r)) + Xi2($/a(r)$/i2(r) 
+ $fi(r)$t2(r)) - F(r)($;ii(r) + ^h2{r) 

where / 'Dra^i{s) = J Pr^ i(s)7'[rQ_i, 0, 1] is a weighted 
functional integral over all the possible configurations of 
the ath chain of lipids or solvent, and V[ra,i,0, 1] is the 
probability distributions of molecular conformations for 
the ath chain, x is the Flory-Huggins interaction pa- 
rameter between the hydrophilic components ( lipid head- 
groups or solvents) and the hydrophobic lipid tails, and 
Xi2 is the head-head or tail-tail pair interactions of dif- 
ferent kinds of lipids. H{r) is the surface field of the 
substrate, and we set[3l| 

{oo, if ^ z < S{x), 

A{i+cos[^(.-s(.))A]} ^ S(x)^z<Six) + e, 
0, if z ^ S{x) + e, 

(7) 

where A describes the field strength at the substrate sur- 
face. 

The mean-field free energy of the system is given 
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il-<Pl-cb2)V' 

dr[xN{tphiir) + tph2ir) + (Psir)){(ptiir) 

+ (pt2{r)) + Xi2N{iphi{r)(ph2{r) + (pti{r)(pt2{r)) 

- H{r)N{(phi{r) + iph2{r) + ifisir) - ipti{r) 

- (pt2{r)) - Whi{r)phi{r) - wti{r)(ptiir) 

- Wh2{r)ph2{r) - Wt2{r)ft2ir) - Ws{r)(fs{r) 
" C(r)(*o(r) - (pM(r) - ipti{r) - iph2(r) 

- (^f2(r) " (Ps{r))], 

where fcs is Boltzmann's constant, T is the tempera- 
ture, and V = J $o(r)cJr is the effective volume of the 
system.dH 0i = ni{vhi+Np^'^)/V and fi = Np^'^/{vhi+ 
Npq^) are the overall volume fraction and the tail vol- 
ume fraction of lipid species A (i = 1) or B (i = 2) 
in the system, respectively, fli, fl2, and fls are the 
single- molecule partition functions of lipid A, lipid B, 
and solvent, respectively, and they are given by fli = 
J 'Dri{s)exp[-Whiiri{l))- dswti{r,{s))] {i = 1,2) and 



r2s = / r>rs(s)exp[-/g dswsirs{s))]. The variables iphi, 
ipti, and ips are the concentration functions of i-type 
headgroups, i-type tails, and solvents, while Whi, wu, and 
Ws are the fields acting on the headgroups, tail segments, 
and solvent segments, respectively, ^(r) is a Lagrange- 
multiplier field which ensures the incompressibility of the 
system. The fields and densities are then obtained by 
minimizing the free energy in Eq.(8), and given as 



wti{r) -■ 

Wh2{r) 

wt2{r) -■ 

Ws{r) = 
Phi{r) = 

Vhi{r) 
V'ti(i-) = 



+Xi2Nji(ph2{r) - jiH{r)N + ji^{r), (9) 
xN{iphiir) + iph2{r) + ipsir)) 
+Xi2Nvt2{r)+H{r)N + ar), (10) 



= X^72(<^ti(r) +Pt2{r)) 

+Xi2Nj2Phi{r) - l2H{v)N + 72^(0(11) 
= X^i^Phii-r) + Vh2{-r) + <Ps(r)) 

+Xi2Npti{v)+H{v)N + ^{v), (12) 
xN{^ti{v) + Vt2{v)) - H{v)N + ^(r), (13) 
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liPhi{v), 

ill dwti{r) 
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(14) 
(15) 

(16) 



Ph2{r) 



2hV 8^2 



^2 dwh2ir) 
= — 92(r,l)g:^(r,l), 

il2 

Ph2{r) = 72Ph2(r), 

<p2f2V 



Vt2(r) 



dsq2{r,s)ql{r,s), 



^2 dwt2{v) 
M2V 

^^2 io 

(1-01-02)1/ dn, 

dws{r) 

(l-(j}i-(j}2)V 



(17) 
(18) 

(19) 



dsqsir, s)qs{r, 1 - sf^O) 



^(^,,.,(r) + ifiuir)) + ip,{r) = $o(r), 



(21) 



where qi{r,s) and ql{r,s){i = 1,2) are the end-segment 
distribution functions of the i-type lipids, and are defined 
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by, [24 



gi(r, s) = j 'DYa,{P[Ya,i,Q,s\5{Y - Ya,i{s)) 

exp{- / dtwu{Yc,.^{t))}, (22) 
Jo 

q|(r, S) = J 'DYa,i'P[Ya^t,S,l]S{Y - Ya^t{s)) 

exp{-u)/ii(rQ,i(l)) - / dtotj(ra^i(t))}. (23) 



(24) 
(25) 



They satisfy the modified difi'usion equations, 
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where a is the statistical segment length. The initial 
conditions are (7,;(r,0) = 1 and ql{Y,l) = exp(— ti;/jj(r)). 
<;s(r, s) has the analogous definition for the solvent, which 
satisfies Eq. (24) with the initial condition gs(r, 0) = 1. 
The resulting self-consistent equations (9) - (21) together 
with equations(22)-(25) can be numerically solved by the 
combinatorial screening method proposed by Drolet and 
Fredrickson. [2^ 

We choose the tail volume fractions in lipid species A 
and B as /i = 0.6 and /2 = 0.45, and the other param- 
eters are fixed to be = 20, N = 50, a = 1, and 
= 1. Thus, we have 71 = 0.67, 72 = 1.22, Vhi — 33, 
and Vii2 = 61, and with these parameters, lipid species 
A and B form a lamellar and a hexagonally cylindrical 
phases in bulk, respectively. [1^ This means that lipid 
species A has a symmetric head/tail shape favoring the 
packing of fiat interfaces, while lipid species B is asym- 
metric, which prefers to curve the interfaces to the lipid 
tails. As lipid species A and B are mixed to put onto a hy- 
drophilic substrate, a mixed lipid bilayer may be formed. 
Here we take the overall volume fractions of lipid species 
A and B as 0i = (/)2 = 0.12, and assume that the repulsive 
head-head/tail-tail interactions between the two kinds of 
lipids are weak(xi2A^ = 3). The strength of surface field 
is taken to be AA^ =1. To facilitate the formation of 
a bilayer structure, we adopt the scenario proposed in 
Ref. 29 by introducing an external field favoring the lipid 
tails in the initial stage of the self-consistent calculation, 
which enforces the bilayer membrane formed along the 
substrate surface profile. When the bilayer is formed, the 
external field is turned off and the calculation continues 
until the equilibrium state of system is reached. 



Results and Discussion 

We now examine the formation of all the possible lipid 
structures by systematically varying the groove depth 
(Hq) and the substrate period (Lq), and the results are 



summarized in Fig. 2. Figure 3 shows a morphology di- 
agram, and a rich variety of lipid domains can be tuned 
by varying Lq and Hq. Here, the height of system is fixed 
to be 40 which is large enough to not influence the mor- 
phology of lipid bilayer, and the lateral length ranging 
from 100 ^ 500 is adjusted to search the stable phases 
by the total free energy minimization. The substrate pe- 
riod changes from Lq = 10 to = 100 and the groove 
depth from Hq = to Hq — 10, which are either compa- 
rable to or several times larger than the natural size of the 
lipids(i?o = cl^/N ~ 7), and therefore, the evolved scale 
of lateral organization is within the nanometer size(i.e., 
from tens to hundreds of nanometers). 

When Hq is vanishingly small, the perturbation from 
the substrate geometry is neglectable, and lipids are 
mixed uniformly on the substrate due to mixing entropy. 
Thus a homogeneously mixed lipid bilayer (HMLB) state 
is formed(Fig. 2(a)). With the increase of Hq, the sub- 
strate becomes slightly rough, and thus lipid species B 
will first occupy the grooves to gain more conformational 
entropy. This is because that, although the substrate is 
relatively fiat, there still exist some perturbations to the 
lipid bilayer from the substrate where the grooves drive 
the formation of lipid B-rich islands (nucleation). As lipid 
B-rich islands initiates from the grooves, a lateral expan- 
sion of the islands on the substrate appears, and even- 
tually leads to the formation of macro-segregated lipid 
bilayer(MSLB) under the help of the lipid A - B in- 
terfacial energy (Fig. 2(b)). When the groove depth 
Hq is increased to modest values where the curvature 
of grooves becomes comparable to the bulk curvature 
of lipid species B, both lipid species B and A can get 
more entropy as they stay in the grooves and on the flat 
steps, respectively. Thus, a component-modulated lipid 
bilayer(CMLB) structure appears in this region, as shown 
in Fig. 2(c). 

When Hq is enough large, the groove curvature may 
be larger than the spontaneous curvature of lipid species 
B at small or modest Lq region, and the lipid distribu- 
tion may become more complicated. In this case, if the 
bilayer structure is still maintained, the bilayer would be 
curved strongly in the regime of grooves, and thus lipid 
species B in the upper leaflet must lose plenty of con- 
formational entropy because the curvature of the upper 
leaflet is in contrast to their bulk curvature. However, 
lipid species B in the lower leaflet prefers to flU in the 
grooves, compared to the case if lipid species A fills in 
the grooves. As a result, the lipid bilayer may be rup- 
tured by the extremely rough substrate: the component- 
modulated lipid domains(CMLD) appears at modest Lq, 
and the coalesced long lipid droplets(CLLD) with spread- 
ing lipid-A droplets is formed at smaller Lq (see Fig. 3). 
The formation of the CMLD and CLLD structures in 
Fig. 2 can be understood by the help of the morpholo- 
gies of thin films on chemically patterned substrates. Ac- 
tually, in the CMLD and CLLD structures, the lipid 
species B almost fill up the grooves, and therefore for 
lipid species A, the substrate now becomes a chemically 




FIG. 2: Summary of all observed morphologies for headgroup of lipid spcies A (iphi) and of lipid species B {Lph2)'- (a) homoge- 
neously mixed lipid bilayer(HMLB), (b) macro-segregated lipid bilayer(MSLB), (c) component-modulated lipid bilayer(CMLB), 
(d) component-modulated lipid domain(CMLD), and (e) coalesced long lipid droplet(CLLD) structures. 




FIG. 3: Morphology diagram for the mixed lipid bilayer sup- 
ported on a geometrically patterned substrate. 



patterned one which is composed of alternating more- 
wettable hard- wall parts and less-wettablc lipid-B parts. 
Thus, the problem can be highlighted by templating of 
thin films on chemically patterned substrates studied by 



Kargupta and Sharma|32|. They examined the effects of 
the periodicity of substrate pattern and the width of the 
more-wettable stripes on the morphology of thin films, 
and found that when the patterned substrate periodicity 
is beyond a critical length, the thin film pattern closely 
replicates the substrate surface pattern, namely dewet- 
ting occurs on every less-wettable parts of the patterned 
substrate, just as shown in Fig. 2(d) where lipid species A 
is almost localized on more-wettable hard-surface parts. 
On the contrary, when the periodicity is smaller than the 
critical length, the dewetting effect may partially be sup- 
pressed, and the lone droplet of thin films spans across 
the remaining less-wettable sites, exactly like the mor- 
phology of lipid-A droplets shown in Fig. 2(e). In real 
biological systems, it was reported [11] that large quan- 
tities of membrane fusion pore sites which are concen- 
trated in cone-shaped nonlamellar lipids such as phos- 
phatidylethanolamine(PE), are simultaneously formed at 
high-curvature regions. This provides further evidence 
implying the validity of these unexpected structures in 
Figs. 2(d) and 2(e). Therefore, the lipid heterogene- 
ity by forming CLLD and CMLD structures may high- 
light the understanding of the relation between the mem- 
brane curvature and fusion, fission, or pore formation of 
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FIG. 4: Free energy contributions from the entropy of lipid 
species A (o), the entropy of lipid species B (•), and interfacial 
energy between the two kinds of lipids(A). (a) Lo = 40; (b) 
Ho = 6. 



membranes. In the large-Z/o region, the curvature of the 
groove become comparable to the spontaneous curvature 
of lipid species B again, and a CMLB phase is presented 
there. 

Importantly, by changing the position of grooves for 
given Hq and Lq, we find that at arbitrary spatial distri- 
bution of grooves, the above conclusions keep unchanged, 
independent of distributed positions of grooves. Actu- 
ally, such a substrate geometry mediated phase-separated 
behavior of lipids provides a controllable mechanism for 
spatial sorting of membrane components by varying the 
roughness of substrates. Furthermore, this is also help- 
ful to understanding how cellular membrane geometry to 
govern the lipid ordering and domain rafts in real biolog- 
ical systems. 

To clarify the underlying mechanism of formed phase 
behavior of mixed lipids, we examine the competition of 
the interfacial energy, conformational entropy, and mix- 
ing entropy of lipids. Figure 4(a) shows the variation of 
free energy F of the system versus Hq by fixing Lq =40. 
When Ho — (the substrate is completely smooth), the 
mixing entropy dominates the lateral organization of the 
two kinds of lipids, leading to the uniform dispersion of 
lipids. When Hq is low (e.g., Hq = 1^^ 3), the inter- 
facial energy dominates the distribution of lipids and a 
MSLB structure appears with lateral growth of lipid B- 
islands located at grooves. For middle values of Hq (e.g., 
Hq — 4~ 7), a CMLB structure is formed by releasing 
the conformational entropy of lipid species B, while at 
the same time, the interfacial energy is increased. Such 
a CMLB structure is mainly contributed to the entropy 
effect of lipid species B. When Hq is large enough (e.g., 
Hq = 8 ~ 10), the formation of CMLD structure where 
the two kinds of lipids separate strongly, largely decreases 
the interfacial energy, while also leads to the loss of mix- 
ing entropy. In Fig. 4(b), we show the variation of free 
energy versus Lq by fixing Hq = 6. The mixed lipid bi- 
layer is either component-modulated structure for small 
Lq as a result of the release of conformational entropy of 



lipid species B or a macro-separated bilayer for large Lq 
where the interfacial energy is greatly decreased. 

It is vitally useful to make a comparison between the 
present work and experimental works. [Ill IHi [3 [3 
First, the appearance of CMLB structure accurately con- 
firms the experimental results observed in Refs. 11 and 
12. We also find that such a CMLB structure formed in 
the modest curvature regions exactly corresponds to the 
nanorafts emerged in the nanocorrugated regions in Ref. 
13 and the spatial modulated structure of alternating the 
lo domains in low curvature and Id domains in high curva- 
ture regions observed in Ref. 14, while the MSLB struc- 
ture presented in the lower curvature regions is just as the 
case that the macroscopic lo domains appear in nanos- 
mooth regions in Ref. 13. Further, our phase diagram 
describing all the observed morphology regimes, which 
greatly enriched the lateral organization of membrane 
component on a geometrically patterned substrate, will 
be very helpful to experimentally probe the lipid distri- 
bution on a wide range of length scales of nanometers. [31 
Finally, by comparing the physical models of the exper- 
imental works and our theoretical work, we easily find 
that although our system is not completely identical to 
those studied by Yoon et al.[ll] and Parthasarathy et 
al.[lj, the physics of these models is consistent, namely 
the curved regions offer some barriers for lipid-A domains 
in the present work and lo domains in Refs. 13 and 14. 
This is because the sphingolipid/cholesterol-rich lo do- 
main and cholesterol-poor Id domain can reasonably be 
represented by symmetric lipid A and asymmetric lipid 
B in our model, respectively. The only difference during 
the lipid organization is that the free energy barrier in 
the experimental works is provided by the bending en- 
ergy of lipid domains, while in our system, comes from 
the loss of conformational entropy of lipids. 



Conclusions 

We conducted a thorough investigation in the lateral 
organization of a two-component lipid bilayer on a ge- 
ometrically patterned substrate by using self-consistent 
field theory. Because of their different molecular shapes, 
the two kinds of lipids prefer to stay in the grooves and 
fiat steps of substrate, respectively. By adjusting the 
groove depth and the substrate period, we obtained a rich 
variety of laterally organized lipid structures. We also 
examine the microscopic mechanism of self-assembly of 
lipids in the mixed bilayer supported on the geometrically 
patterned substrate, and find that the resulting phase be- 
havior of mixed lipids is dominated by the competition of 
the interfacial energy, conformational entropy, and mix- 
ing entropy of lipids. The present study confirmed recent 
experimental results in Refs. 11 and 12, and qualitatively 
agreed with the experimental findings of Refs. 13 and 14. 
Most importantly, we also predict, for the first time, the 
formation of the CLLD and CMLD structures in large 
curvature regions, which may provide theoretical insight 
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into future experiment with the improvement of visual- No. 20674037, and No. 10574061. 
ization techniques towards smaUer scales of geometry and 
help to understand the microscopic mechanism of fusion 
pore formation in biomembranes. 
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